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ABSTRACT 

Following two years of complete occultation of both stars by its opaque circumbinary ring, the 
binary T Tauri star KH 15D has abruptly brightened again during apastron phases, reaching 7 = 15 
mag. Here, we show that the brightening is accompanied by a change in spectral class from K6/K7 
(the spectral class of star A) to ^Kl, and a bluing of the system in V — I by about 0.3 mag. A radial 
velocity measurement confirms that, at apastron, we are now seeing direct light from star B, which is 
more luminous and of earlier spectral class than star A. Evidently, the trailing edge of the occulting 
screen has just become tangent to one anse of star B's projected orbit. This confirms a prediction of 
the precession models, supports the view that the tilted ring is self-gravitating, and ushers in a new 
era of the system's evolution that should be accompanied by the same kind of dramatic phenomena 
observed from 1995-2009. It also promotes KH 15D from a single-lined to a double-lined eclipsing 
binary, greatly enhancing its value for testing pre-main sequence models. The results of our study 
strengthen the case for truncation of the outer ring at around 4 AU by a sub-stellar object such as an 
extremely young giant planet. The system is currently at an optimal configuration for detecting the 
putative planet and we urge expedient follow-up observations. 
Subject headings: circumstellar matter — protoplanetary disks — binaries: close 



1. INTRODUCTION 

KH 15D is a young binary system composed of similar, 
but not identical, low-mass pre-main sequence stars in 
an orbit of eccentricit y ~0.6 with a period of 48.37 days 
([Johnson et al.l 12004 IWinn et al.l 12004 Hamilton etaTJ 
120051 ) . The binary orbit is viewed nearly edge-on 
and is embedded in an accretio n disk from which a 
well-collimated out f low emerges ([Hamilton et al.l 1200 



Deming et all 120041: iTokunaga et al.l l2004t iMundt et al.l 



20101 The estimated age of the stars is 3 x 10' 



based on their location in the color-magnitude diagram 
of NGC 2264, and the total system mass is ~1.3 M© 
([Hamilton et al.ll2001l : [Johnson et al.ll2004 ). A thin cir- 
cumbinary ri ng of evolved solids has precipitated from 
the gas disk ([Lawler et al.l |2010|) within the terr estrial 
zone (1-4 AU) of this system (|Herbst et al.ll200cl . The 
ring reveals itself as an op aque screen w ith a razor-sharp 
edge (< 0.1 stellar radii; iHerbst et al.l 2010) at optical 
and near-infrared wavelengths. For decades, the lead- 
ing edge of this structure has been slowly moving across 
the binary orbit, apparent ly resulting from precession of 
the slightly inclined rin g ([Chiang fe Murray- Clavl 120041 : 
iWinn et al.l 12004 I2006D . One star's orbit (designated 
star B) was com pletely occult ed in 1995 and the other 
(star A) in 2009 (|Herbst et afll2010l ). 
The tidal influence of a binary on a misaligned disk 



is known to cause rigid-body precession given that 
the disk can communicate at length-scales comparable 
to the inner disk radius and time-scales close t o the 
local precession perio d (lArtvmowicz fc Lubowl 11994 
Papaloizou fc Terqueml 119951: ILarwood fe Papaloizoul 



1997|) . The variable torque applied by a non-coplanar 



binary will truncate the inner disk radius and induce 
a gradient of inclinations across concentric annuli in 
the disk, producing a warp. The specific geometry of 
the warp, such as whether the gradient increases or 
decreases with radius, depends upon the mechanism 
responsible for enforcing an alignmen t of t he nodes 
of the annuli. iChiang fc Murrav-Clavi ([2004 ) propose 
two means by which the KH 15D circumbinary ring 
maintains such an alignment - either by self-gravity 
of the ring particles or by thermal pressure normal to 
the ring plane - with a unique warp geometry resulting 
in each case. In such warped, precessing disks, if the 
outer radius is allowed to become arbitrarily large, 
then the disk precession period will approach infinity. 
Thus, for nodal precession to occur, the outer disk 
must also be confined in radial extent (i.e. it forms a 
small, separate ring), s uch as by a shepherding body 
(IGoldreich fc Tremainel [19821 IChiang fc Murrav-Clavi 
12004) or a sufficiently rapid d ecrease in surface density 
(Larwood & Papaloizou 1997). The misalignment of the 
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binary orbital plane with the KH 15D ring is believed 
to induce precession in its ring and, further, the system 
exhibits distinctive eclipses due to this non-coplanar 
arrangement. 

Although identified as a va riable star in the late 1960s 
(jBadalian fc Erastoval I1970T ). i nterest in KH 15D be- 
gan with the discovery in 1995 (|Kearns fc Herbstl [19981) 
that the star essentially winked on an d off by several 
magn itudes on a period of 48.37 days ([Hamilto n et al.l 
2005), now known to be the orbital period of the bi- 
nary (jJohnson et al.ll2~004l ) . These extreme brightness ex- 
cursions started then because precession of the ring had 
moved one of its razor-sharp edges fully across the orbit 
of one member of the binary (star B), leaving the other 
(star A) to periodically rise and set with respect to that 
horizon as seen from Earth. From 1995 to 2009, the wink- 
ing continued, with the system being increasingly more 
"off" as star A spent more and more of each orbital cy- 
cle be low the horizon ([Hamilton et alf2005t iHerbst et al.l 
120101) . In 2010 and 2011, the system was fainter at all 
times because neither star appeared above either edge. 

Both stars are believed to be less massive than the 
Sun, with star A haying a spectral type of K6/K7 
([Hamilton et al.ll2001b lAgol et al.ll2004f ). The prediction 
that star B is slightly more massive and of earlier type 
has heretofore only been inferred from the fact that it 
is a brighter star, as evidenced by archival photometry 
(iWinn et al.ll2003t [Johnson fc Wiimll2004l Ijohnson et all 
120051 ) and bv its one appearance in the modern photo- 
metric record (see Fig. [IJ prior to 2012. 

2. OBSERVATIONS AND DATA ANALYSIS 

We obtained images of KH 15D on over three hun- 
dred nights between October 2010 and March 2012, using 
the ANDICAM instrument on the 1.3m telescope oper- 
ated by the SMARTS consortium at Cerro Tololo Inter- 
American Observatory (CTIO). The /—band photome- 
try is in addition to an observing record from multiple 
observatories dating back to 1995 and the V — I color 
information has been on record for nearly a decade. The 
optical images consisted of four exposures, each lasting 
150 seconds. We used IRAF scripts to flat field, sky sub- 
tract and combine dithered images, and the package phot 
to complete the aperture photometry. The median error 
is 0.01 mag for / < 17.5 mag and increases to 0.05 mag 
at / = 18.5 mag. 

We obtained cross-dispersed (XD) spectroscopy with 
the Gemini Near Infrared Spectrograph (GNIRS) on UT 
January 5, 2011, when the stars were both near apas- 
tron but fully obscured, with / = 16.95 mag. The XD 
data were taken in the GNIRS short focal length camera 
mode with wavelength coverage from 0.9 to 2.5 /im, a 
slit width of 0.3", and a 31.7 line mm -1 grating (1600 < 
R < 1700). The processing used the gemini.gnirs IRAF 
package Beta Version 1.11. The data products consist of 
a set of calibration images, a spectrum of KH 15D and 
a spectrum of HD 59374, an F8 V star observed for the 
purposes of characterizing the sensitivity function. The 
integration times were less than four seconds for all of 
the calibration and standard exposures. The observa- 
tions of the target itself include 20 frames, each of 300- 
second duration. The spectra were flat-fielded, assigned 
a wavelength solution by comparison with an argon arc- 
lamp spectrum and corrected for spatial distortions and 
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Fig. 1. — The long-term light curve of KH 15D showing its dra- 
matic re-brightening in 2012. The one measurement in 1995 show- 
ing the system at / = 14 mag was obtained near periastron, when 
star B made its last appearance above the leading edge of the 
screen. 
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Fig. 2.— Light curves from 2010 and 2011 (x's) and 2012 (solid 
circles) phased with the orbital period of 48.37 days. The dramatic 
brightening that occurred around apastron (phase = 0.5) in 2012 is 
apparent and represents the appearance of star B partially above 
the ring edge. 

detector artifacts. The details o f the standard GNI RS 
reduction procedures are given in IWinee et al.l ((20091) . 

3. RESULTS 

3.1. Lightcurve 

The /-band light curve of KH 15D from 1995 to 
present is shown in Fig. [T] Each vertical strip represents 
an annual observing season, culminating in October- 
March 2012. The full on-off eclipse behavior, caused 
by star A disappearing behind the leading edge of the 
occulting screen during each orbit, was present up until 
about 2006. The maximum system brightness then grad- 
ually diminished, because star A never fully rose above 
the screen edge. The system spent increasingly more 
time in complete obscuration until 2010, when both stars 
were fully occulted at all phases and the system never ap- 
peared brighter than / = 16.3 mag. Interestingly, even 
when fully occulted, the system varied by more than two 
magnitudes on the orbital period, being brighter near 
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Fig. 3. — A comparison of the colors in 2012 to four years prior. The set of points which collect around a relatively blue median color of 
V — I = 1.25 mag represent star B and those that collect around a redder median color of V — I = 1.56 mag represent star A. 



apastron and fainter near periastron, due to scattered 
light (see Fig.©. 

The re-brightening of the system that occurred in 2012, 
reaching / = 15 mag near apastron (phase = 0.5), is ap- 
parent on Fig. [T]and highlighted on Fig.© where we show 
light curves for each of the past three seasons, phased to 
the orbital period. A "star rise" is evident during one 
cycle in 2012 at a phase of 0.35 and a "star set" dur- 
ing a different cycle is seen at phase 0.65. The maximum 
brightness near / = 15 mag is less than the known bright- 
ness of star A (I — 14.5 mag) or star B (J = 14.0 mag) 
and, therefore, most likely represents a star still partially 
occulted by the ring. This is analogous to the lightcurve 
in the 2008 season, when star A was partially occulted. 
We now show that this is not star A reappearing, but 
star B appearing from behind the trailing edge of the 
occulting screen. 

Fig.[3]shows the V — I color as a function of brightness, 
comparing the photometry from 2008 and 2012. The 
horizontal lines in the figure show the median color for 
each season when the system is bright, i.e., after star rise. 
In 2008, we found V — I = 1.56 mag, characteristic of 
the K6/K7 star (A). In 2012, we measure V — I = 1.25 
mag, significantly different and characteristic of a hotter 
star. Correcting for foregr ound reddening [E{B — V) — 
0.15 {Hamilton ct al. 200l|), corresponding to E(V — I)= 
0.21], we find the intrinsic color of the (partly) visible star 
in 201 2 to be (V — I)n = 1-04 mag, characteristic of a Kl 
giant (Bcssell 1979). Note that there is no evidence of any 
reddening accompanying the large extinction suffered by 
the stars, either in 2012, or in any year between 1995 and 
2009, when star A was undergoing regular rise and set 
events. 

3.2. Spectral Classification 

The Gcmini/GNIRS spectrum obtained in January 
2011, when the system had I — 16.95 mag, supports our 
photometric result that a bluer star, presumably star B, 
recently became a significant contributor to the light of 
KH 15D. Fig. [4] shows the spectrum spanning five of the 
XD orders (0.9 to 2.5 fim). Obvious features are labeled; 



the molecular hydrogen emission featur es most likely 
arise from the shocked amb ient medium (|Deming et al.l 
I2004t iKusakabe et~aTl I2005D. the He! line from the fun- 
nel accretion flow (jFischer et al.ll20(M K and the CO band 
heads from the stellar photosphere itself. The spectrum 
was flux calibrated by integrating over the 2MASS fil- 
ter profiles and scaling the observations to fit. The blue 
end was calibrated using the overlap region with the J 
band. The source variability introduces error into the 
flux normalization, but does not affect the spectrum's 
shape nor relative line strengths, which we now use to 
classify the source by visually comparing i t to template 
spect ra (|Pickleslll99"8l : iCushing et "al1l2005t iRavner et all 
I2009D . In Fig.H KH 15D is in black and the red (top) 
and blue (bottom) spectral templates correspond to HD 
25975, a Kl III star, and to HD 237903, a K7 V star, 
respectively. The continuum slope is inconsistent with a 
spectral class as late as K6/7, as expected of star A: the 
KH 15D and K7 V star continua cease to overlap blue- 
ward of 1.3 firn. Matching the continuum slope of the 
full wavelength region in question requires a contribution 
from a star of type Kl±0.5. The set of lines between 1.1 
and 1.2 /jm and also those just before 2 /xm would be 
deeper, were the K7 star the only source. This spectrum 
was obtained near apastron but at a time when both 
stars were fully obscured so it represents the combined 
scattered light. We note that since th e stars interact 
with one another and their disk (see e.g. lHamilton et al.1 
2012), they might defy a clear luminosity classification, 
but as they are still contracting to the main sequence, 
one expects them to effectively be sub-giants or giants, 
in agreement with our findings. 

3.3. Radial Velocity 

To confirm our interpretation that it is star B that 
has caused the system to brighten, we obtained a 10- 
minute exposure of KH 15D with the Keck High Reso- 
lution Echelle Spectrometer (HIRES) on UT March 13, 
2012, when the system was near its maximum bright- 
ness. The radial velocity (RV) was measured by cross- 
correlation with a K0 main sequence star (HD 185144) 
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Fig. 4. — The GNIRS spectrum of KH 15D (in black) compared to template spectra of type Kl III (top, red) and K7 V (bottom, 
blue). Some of the absorption features unique to the KH 15D spectrum, including the prominent ones at 1.275 and 1.546 fim, are as yet 
unidentified. 



yielding a heliocentric value of +39.9 km s 1 . At this 
orbital phase (0.41), the heliocentric velocity of star A 
is predicted to be -0 .1 km s" 1 (Model 3 , I Winn et"aH 
2006) or +2 km s" 1 (| Johnson et al.ll200l . We caution 
the reader that the systemic velocity of the KH 15D 
system is not yet well determined due to the difficulty 
in obtaining spectra suitable for cross-correlation anal- 
ysis while deep in eclipse. Ho wever, two valu es exist in 
the literature ( +14.7 km s " 1 , Uohnson et al.l 2004 and 
+ 18.6 km s-\|W inn et all 2006), which are consistent 
with th e radial velocities of th e stars in the NGC 2264 
cluster dSoderblom et al.iri999P ). We adopt the value of 
I Winn et al.l (|2006l ). finding that relative to the center of 
mass, star B is moving with a velocity of +21.3 ± 1 km 
s _1 , while the predicted velocity for star A is —18.7 ± 1 
km s _1 . The uncertainty represents the standard devi- 
ation of the mean RV measured from all echelle orders. 
These velocities suggest that the two stars are of similar 
mass, however, that star A, not star B, is the more mas- 
sive of the pair. While this seems contradictory (as the 
brighter, hotter star B should be the more massive one), 
we urge discretion regarding this point since the value 
used for star A is model-dependent. Clearly, once more 
RV measurements are obtained as the orbit of star B is 
revealed, the models can be refined. More importantly, 
the RV measurement confirms that we are now observing 
star B at apastron, not star A. 



4. DISCUSSION 

4.1. Confirmation of the Precession Model 

The precessing circumbinary-ring model has provided 
a success ful framework for interpreting observations of 
KH 15 D (IWinn et al.ll200l 120061: IChia ng fc Murrav-"Clavl 



120041: iSilvia fc Agoll 120081) . A prediction is that the 
eclipse behavior exhibited by star A from the early 1990s 
through 2009 should repeat itself with star B as the vis- 
ible star. The developments reported here confirm that 
prediction and provide strong support for this scenario. 
Fig. [S] is a schematic drawing of the system, as we be- 
lieve it appears in 2012. For the first time, we are able 
to locate the trailing edge of the ring which provides im- 
portant new information to improve the next generation 
o f models. 

I Winn et al.1 (|2004L l2006h employed archival records to 
address the secular evolution of KH 15D by projecting 
the binary orbit onto the plane of the sky and treating 
the precession of the ring as a sharp, semi-infinite edge 
(that is, the location of leading edge is finite and the 
trailing edge extends to infinity) that marches across the 
projected binary orbit over time. The authors predicted 
the eventual re-emergence of star B as a natural result 
of the advancing edge. As the model had no way to 
differentiate between viewing angles with respect to the 
occulting feature nor its thickness, it provided no definite 
prediction regarding whether both binary orbits would 
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ever be occulted at the same time or, if so, when star B 
might reappear. Having the location of only one edge 
of the screen specified also complicated the modeling of 
the scattered light. The discovery of the location of the 
trailing edge should lead to a much improved version of 
this model which will ultimately quantify most of the 
i mportant physical characteristi cs of the system. 

IChiang fc Murray- Clavl (|2004D presented a first-order 
dynamical theory of the KH 15D ring noting that its 
rigid precession can be enforced either by thermal pres- 
sure gradients or by ring self-gravity. They favored self- 
gravity and their case is strengthened by the observations 
contained herein. The future light curve of KH 15D pre- 
dicted by IChiang fe Murrav-Clavi (|2004T> for the thermal 
pressure case (Model 1) differed substantially from that 
for self-gravity (Model 2) because the direction of warp- 
ing is different. In Model 1, the disk is closer to the 
plane at large distances and more inclined closer in and 
the predicted duration of full occultation is quite long 
(many decades). In Model 2, the disk is more inclined 
at larger distances and the duration of full occultation 
of both orbits is short or non-existent. We now know 
that the observed light curve had a short duration (^2 
years) of full occultation, qualitatively in much better 
agreement with the self-gravity model, although quanti- 
tatively in agreement with neither. 

While recent data will enable properties such as the 
ring geometry, orientation and scattering properties to 
be refined, the mounting support for rigid nodal preces- 
sion underscores that the occulting ring of solids must be 
confined at its outer borders. Therefore, the inferred ra- 
dial dimensions of the ring should not change drastically 
from past estimates. 

4.2. Remaining Questions and Future Observations 

The question r emains as to what has truncated the ring 
at around 4 AU. IChiang fc Murrav-Clavl ()2004D suggest 
a planet as the most likely cause. Being at or near the 
ice line in the system it is possibly a giant planet in for- 
mation. Color excursions to values redder than star A 
have long been known in KH 15D during minimum light 
and may find their interpretation in the light of this pu- 
tative third body. Such an extremely young giant planet, 
if sufficiently massive, might be self-luminous enough to 
be detectable in the near or mid-infrared during phases 
when both stars are occulted. It should be sought im- 
mediately, as the minimum brightness of the system is 
already beginning to gradually increase due to the ad- 
vancement of the trailing edge (see Fig. [1]). 

Fig. [U shows that the abrupt change in slope that oc- 
curs during star rise and star set continues, suggesting 
that the trailing edge of the ring will be as sharp as the 
leadi ng edge: namely, m uch smaller than a stellar ra- 
dius (jHerbst et al.l l2010h . Resultantly, the eclipses that 
will ensue will allow us to determine the radius of star B 
with the same precision that can be achieved in a typ- 
ical eclipsing binary. We can now obtain the RV curve 
of star B, and so will constrain the masses, radii, lumi- 
nosities and effective temperatures of the stars. KH 15D 
will soon be the first known double-lined eclipsing binary 
young enough to still be embedded in an accretion disk, 
making it of paramount importance for testing models of 
pre-main sequence evolution. 

As Fig. [5] shows, the physical extent of the obscuring 
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Fig. 5. — A schematic representation of the KH 15D system as 
projected on the sky. 



portion of the ring is slightly wider than the projected 
size of the binary orbits. What sets this width remains 
unclear to us. Multiple possibilities exist but it is well 
beyond the purview of this observational paper to as- 
sess them here. The KH 15D spectral features of non- 
stellar origin will be addressed in subsequent work, and 
additional spectra should be obtained to verify these fea- 
tures and probe more sight lines in the ring. Over the 
coming decade we will have the opportunity to use the 
trailing edge for natural coronography to study the mag- 
netosphere and jet-launching zone of star B in the same 
detail as was done for star A. 

5. SUMMARY 

The reappearance of star B provides a dramatic confir- 
mation of the precessing-ring model, demonstrating that 
such structures must be a common feature of circumbi- 
nary and, presumably, circumstellar transition disks. It 
opens the door for a more complete study of the ring 
and what its role in planet formation may be. This 
is now more pertinent than ever, as the Kepler mis- 
sion c ontinues to discover circumbinary pl anetary sys- 
tems (|Dovle et al.l [20111: IWelsh et al.l [2012]) . Follow-up 
observations of KH 15D should be conducted forthwith to 
uncover the putative planet confining the circumbinary 
ring. The emergence of star B has transformed KH 15D 
from a single-lined spectroscopic binary to a double-lined 
eclipsing system, with all that that entails for precise 
measurement of fundamental astrophysical properties. 



Acknowledgments 

We thank Geoff Marcy for helping to arrange the 
Keck/HIRES observations. We are grateful to the ob- 
servers and personnel at the SMARTS and GNIRS fa- 
cilities. This work was partially supported by a NASA 
grant to W.H. through the Origins of Solar Systems pro- 
gram. Based on observations obtained at the Gemini 
Observatory via queue program GN-2010B-Q-52. Gem- 
ini is operated by the Association of Universities for Re- 
search in Astronomy, Inc., under an agreement with the 
NSF on behalf of the Gemini partnership: the National 
Science Foundation (U.S.), the Science and Technology 
Facilities Council (U.K.), the National Research Council 
(Canada), CONICYT (Chile), the Australian Research 
Council (Australia), Ministerio da Ciencia, Tecnologia e 
Inovagao (Brazil) and Ministerio de Ciencia, and Tec- 
nologia e Innovation Productiva (Argentina) . 



6 



Capelo et al. 



REFERENCES 



Agol, E., Barth, A. J., Wolf, S., & Charbonneau, D. 2004, ApJ, 
600, 781 

Artymowicz, P., & Lubow, S. H. 1994, ApJ, 421, 651 
Badalian, H. S., & Erastova, L. K. 1970, Astronomichcskij 

Tsirkulyar, 591, 4 
Bessell, M. S. 1979, PASP, 91, 589 

Chiang, E. I., & Murray-Clay, R. A. 2004, ApJ, 607, 913 
Cushing, M. C, Rayner, J. T., & Vacca, W. D. 2005, ApJ, 623, 
1115 

Deming, D., Charbonneau, D., & Harrington, J. 2004, ApJ, 601, 
L87 

Doyle, L. R., Carter, J. A., Fabrycky, D. C. et al. 2011, Science, 
333, 1602 

Fischer, W., Kwan, J., Edwards, S., & Hillenbrand, L. 2008, ApJ, 
687, 1117 

Goldreich, P., & Tremaine, S. 1982, ARA&A, 20, 249 
Hamilton, C. M., Herbst, W., Mundt, R., Bailer-Jones, C. A. L., 

& Johns-Krull, C. M. 2003, ApJ, 591, L45 
Hamilton, C. M., Herbst, W., Shih, C, & Ferro, A. J. 2001, ApJ, 

554, L201 

Hamilton, C. M., Herbst, W., Vrba, F. J. et al. 2005, AJ, 130, 
1896 

Hamilton, C. M., Johns-Krull, C. M., Mundt, R., Herbst, W., & 

Winn, J. N. 2012, ApJ, 751, 147 
Herbst, W., Hamilton, C. M., Leduc, K. ct al. 2008, Nature, 452, 

194 

Herbst, W., LeDuc, K., Hamilton, C. M. et al. 2010, AJ, 140, 2025 
Johnson, J. A., Marcy, G. W., Hamilton, C. M., Herbst, W., & 
Johns-Krull, C. M. 2004, AJ, 128, 1265 



Johnson, J. A., & Winn, J. N. 2004, AJ, 127, 2344 
Johnson, J. A., Winn, J. N., Rampazzi, F. ct al. 2005, AJ, 129, 
1978 

Kearns, K. E., & Herbst, W. 1998, AJ, 116, 261 
Kusakabe, N., Tamura, M., Nakajima, Y. ct al. 2005, ApJ, 632, 
L139 

Larwood, J. D., & Papaloizou, J. C. B. 1997, MNRAS, 285, 288 
Lawler, S. M., Herbst, W., Redfield, S. et al. 2010, ApJ, 711, 1297 
Mundt, R., Hamilton, C. M., Herbst, W., Johns-Krull, C. M., & 

Winn, J. N. 2010, ApJ, 708, L5 
Papaloizou, J. C. B., & Terquem, C. 1995, MNRAS, 274, 987 
Pickles, A. J. 1998, PASP, 110, 863 

Rayner, J. T., Cushing, M. C, & Vacca, W. D. 2009, ApJS, 185, 
289 

Silvia, D. W., & Agol, E. 2008, ApJ, 681, 1377 

Sodcrblom, D. R., King, J. R., Siess, L., Jones, B. F., & Fischer, 

D. 1999, AJ, 118, 1301 
Tokunaga, A. T., Dahm, S., Gassier, W. et al. 2004, ApJ, 601, 

L91 

Welsh, W. F., Orosz, J. A., Carter, J. A. et al. 2012, Nature, 481, 
475 

Wingc, C, Riffel, R. A., & Storchi-Bergmann, T. 2009, ApJS, 
185, 186 

Winn, J. N., Garnavich, P. M., Stanek, K. Z., & Sasselov, D. D. 

2003, ApJ, 593, L121 
Winn, J. N., Hamilton, C. M., Herbst, W. J. et al. 2006, ApJ, 

644, 510 

Winn, J. N., Holman, M. J., Johnson, J. A., Stanek, K. Z., & 
Garnavich, P. M. 2004, ApJ, 603, L45 



